Controlled vaporization of the superconducting condensate in cuprate 
superconductors sheds light on the pairing boson. 
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We use ultrashort intense laser pulses to study superconducting state vaporization dynamics in 
La2-iSra;Cu04 (x = 0.1 and 0.15) on the femtosecond timescale. We find that the energy density 
required to vaporize the superconducting state is 2 ± 0.8 K/Cu and 2.6 ± 1 K/Cu for x = 0.1 
and 0.15 respectively. This is significantly greater than the condensation energy density, indicating 
that the quasiparticles share a large amount of energy with the boson glue bath on this timescale. 
Considering in detail both spin and lattice energy relaxation pathways which take place on the 
relevant timescale of ~ 10 -12 s, we rule out purely spin-mediated pair-breaking in favor of phonon- 
mediated mechanisms, effectively ruling out spin-mediated pairing in cuprates as a consequence. 

The study of noncquilibrium phenomena in superconductors has been an important topic of condensed matter 
physics since the 1960's and the fact that intense laser pulses can non-thermally destroy the superconducting state 
has been known for a long time pQ . After the discovery of high-temperature superconductivity in cuprates and the 
simultaneous rapid development of ultrashort pulsed lasers, real-time studies of quasiparticle (QP) dynamics became 
possible using pump-probe experiments [3J 0] [5]. Recent developments in phenomenological modeling and new 
systematic experimental studies of the non-equilibrium optical response [6j [7j [8j [9] have enormously improved our 
understanding of the dynamics of photoexcited QPs on short timescales, and the response of the superconducting 
state to weak pulsed laser excitation in cuprates can now be unambiguosly identified on the femtosecond timescale 
[3 13 HI H0l HU H21 ESI EH HSl EH]. Rothwarf and Taylor originally proposed their phenomenological model for QP 
recombination in the framework of phonon-mediated pairing[6, 8 , but spin-fluctuation mediated recombination (which 
is relevant for the cuprates) is not excluded by their model. Both high-energy phonons and spin excitations could in 
principle mediate QP recombination, and the non-equilibrium studies so far did not directly reveal the pairing boson. 

In this work we specifically address the question of the mediating boson by carefully measuring and analyzing the 
energy and time needed to destroy (vaporize) the superconducting condensate in L^-^Sr^CuC^ with x — 0.1 and 
0.15. We compare the measured vaporization energy with the thermodynamically measured condensation energy, and 
find that a substantial amount of energy is temporarily stored by the glue boson bath during the vaporization process. 
Since spin and lattice subsystems have vastly different heat capacities, this places significant constraints on the type 
of bosonic bath which can mediate pairing. By carefully considering the energy relaxation pathways associated with 
pair-breaking dynamics we are able to conclude which bosonic excitations are involved in the destruction of the 
condensate, shedding light on the pairing boson responsible for superconductivity in these materials. 

Outlining the sequence of events in our experiments microscopically, the laser pump pulses first excite electrons 
from occupied to unoccupied states within 1.5 eV of the Fermi level. Immediately afterwards, in an avalanche QP 
multiplication process which is well-studied in metals, as well as cuprates, the photoexcited carriers relax to states 
near the Fermi energy via mfraband electron-electron scattering, occuring on a typical timescale r e _ e < 50 fs |17j . 
and scattering with phonons - preferentially interacting with those phonons which are most strongly coupled to the 
QPs [18 - resulting in significant non-equilibrium QP and phonon populations within ~100 fs of photoexcitation. The 
next relaxation step, QP recombination across a superconducting energy gap (or pseudogap) with the emission of a 
boson with energy > 2A, takes significantly longer, and is typically described very well by the Rothwarf and Taylor 
model 6 , 8 9 . The model does not directly identify the pairing boson, but two crucial parameters in the model do 
depend on the electron-boson interaction, namely the characteristic bare pair-breaking rate r\ and recombination rate 
R which define the pair-breaking and QP recombination timescales respectively. [BJ |7J [5] Importantly, r\ and R can 
be determined from the vaporization time T r .[7] 

The important feature of pump-probe experiments is that the transient density n p of the photoexcited (PE) QPs 
accumulated at the gap edge is probed in real time by a time-delayed excited state absorption process [T3] . The 
transient change in reflectivity AR(t), of the probe pulse (which for small AR is linearly proportional to the transient 
photoinduced absorption [19]) is thus directly proportional to n p (t). This makes it possible to detect when the 
superconducting condensate is destroyed. [5] [T3] 

In the experiments described here, we perturbed the superconducting condensate in L^-xSr^CuC^ with 50 fern- 
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tosecond laser pulses. The experiments were performed on freshly cleaved surfaces of high quality La2_ x Sr a; Cu04 
(x = 0.1 and 0.15) single crystals with T c = 30 K and 38 K respectively [20]. The laser pulses were linearly polarized 
and incident along the c axis of the crystal with a wavelength of A — 810 nm (~1.5 eV). We used a Ti:Sapphirc 
oscillator and a 250 kHz amplifier to cover the range of excitation fluences from T ~ 4 x 10 -2 /iJ/cm 2 to 100/iJ/cm 2 . 
The pump and probe beam diameters were measured accurately with a pinhole and the absorbed energy density 
was accurately determined [21j. The low laser repetition rate of our laser ensured that there was no heat buildup 
between pulses even with the highest fluences used, and the temperature increase due to the laser was found to be 
less than 2 K (which can also be seen from a comparison of the T c measured optically with the T c from susceptibility 
measurements) . 

The photoinduced reflectivity change AR/R as a function of time delay for different T is shown below T c (T — 4.5 
K) in Fig. la) and above T c (T = 32 K) in Fig. lb) for x = 0.1 (the data for x — 0.15 is qualitatively the same). 
Below T c (Fig. la)) we identify two relaxation processes with very different dynamics, which we label as A and B. 
Signal B is present from low T to well above T c (Fig. lb), and disappears gradually above the so called pseudogap 
temperature T* . In agreement with many previous low- J 7 experiments [9j EJ [T2J [13] , it is assigned to the carriers 
recombining across the pseudogap. Signal A is visible strictly only below T c and is - in accordance with previous works 
[111 116]- assigned to QP recombination across the superconducting gap A S (T), and has a relaxation time typically 
ta > lOps at 4.5K [11 . The rise-time r r = 0.8 ± 0.15 ps of the superconducting signal AR/Ra is the time required 
for the QP population to build up [8] by pair-breaking from the condensate. 

Examining Figure 1 in more detail, we see that at low T and T signal A is dominant. As fluence is increased, the 
amplitude of signal A first increases with T and then starts to saturate for T above « 12 /xJ/cm 2 . As signal A starts 
to saturate, signal B starts to become more visible, and above the saturation threshold of signal A, a linear increase 
of the amplitude of signal B with increasing T becomes clearly apparent. 

The maximum amplitudes of AR^/R and ARb/R are shown in Fig. 2a) for x — 0.1 and 0.15 as a function of J- . 
We see that AR/\/R is linear at low fluence for T < 8/iJ/cm 2 . Above 8fi J/cm 2 , the signal amplitude departs from 
linearity, indicating an onset of saturation associated with vaporization of the condensate. AR/\/R soon saturates 
and becomes constant for T > 18/iJ/cm 2 (up to the highest fluences measured). In contrast, ARs/R is linear with 
F both below and above T c . 

To accurately determine the vaporization threshold, we carefully take into account the optical penetration depth 
X op for the pump and the probe beams and their spatial profile [21] . From fits of the measured dependence of AR/\/R 
on T to the function provided by a straightforward model calculation (shown in Fig. 2a)), [21] we obtain values for 
the threshold vaporization fluence at 4.5 K: Tt — 4.2 ± 1.7/iJ/cm 2 for x — 0.1 and Tt — 5.8±2.3/iJ/cm 2 for x — 0.15. 

In Fig. 2b) we plot the T- dependence of AR^/R for several excitation levels for the x = 0.1 sample. As expected, 
for T > Tt, the T-dependence of AR/\/R does not depend on T, since full vaporization is achieved at all T < T c . 
Near the threshold, for T = 7/iJ/cm 2 , only partial vaporization is evident and the amplitude AR&/R merges with the 
high fluence data only as T — > T c . We can understand the T-dependence of the AR^/R by considering the difference 
in reflectivity between the superconducting state and the normal state. The induced change in reflectivity for fluences 
above the vaporization threshold A s — \ ^\ J r >J r T is proportional to er" — af, where er™ and erf are real parts of the 
complex conductivity in the normal and superconducting states, respectively. Using the Mattis-Bardeen formulae [31 
it follows that [21] : 



, , , 2A(T) , flA7huj\ 

where huj is the photon energy and A(T) the T-dependent gap. Using A(T) = Ao ^1 — (T/T c ) 2 ^j (Aq is gap at 
K), which was previously found to describe A(T) in cuprate superconductors [32] . a very good agreement between 



Eq.(20l and the data for T > Tt is obtained (see Figure 2 b). 

Let us now examine the energy relaxation pathways on the pair-breaking timescale of ~ 1 ps. Phonons released 
during this time need at least A op /u s ~ 30 ps to escape from the excited volume, v s being the velocity of sound. 
The characteristic QP diffusion time from the excitation volume is also of the order of ~ 100 ps, calculated using the 
measured QP diffusion constant for very clean samples of YBa2Cu3C>6.5 at 4 K [24] . Therefore we can conclude that 
the absorbed optical pulse energy cannot diffuse or escape, and remains in the excitation volume on the timescale of 
1 picosecond. 

Next, let us analyze the microscopic energy relaxation processes within the excitation volume in more detail. The 
energy densities in the excitation volume at vaporization threshold for x — 0.1 and x — 0.15 shown in Fig. 2 a) are 
U p = T T /(\opkB) = 2.0 ±0.8 K/Cu and 2.6 ± 1.0 K/Cu respectively (using \ op = 150 nm at 810 nm [30])). Both are 
significantly higher than the thermodynamically measured condensation energies extracted from specific heat data, 
which are U c /ks = 0.12 K/Cu for x — 0.1, and U c /ks = 0.3 K/Cu for x — 0.15 [5S]. The ratio of the two energies 
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FIG. 1: (Color online). The photoinduced reflectivity AR/ R in Lai.gSi'o.iCuCU (T c = 30 K) taken at various photoexcitation 
fluences a) below and b) above T c . The data above T c are normalized with respect to T and fully overlap, showing that the 
response is linear. Below T c the two distinct relaxation components are marked as A and B. 
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FIG. 2: (Color online), a) The maximum amplitude AR&/R at 4.5 K for x = 0.1 (empty circles) and x — 0.15 (full circles), and 
ARb/R and at 4.5 K (squares) and 32 K (triangles) for x=0.1 in La2- a: Sr :! ;Cu04. The arrows mark the vaporization thresholds 
Tt = 4.2 ± 1.7/xJ/cm 2 and Tt = 5.8 ± 2.3/xJ/cm 2 for x=0.1 and x=0.15 respectively obtained from the fit [2TJ (lines), b) The 
T-dependence of ARa/R for x=0.1. The dashed line is a fit to the data using Eq.(l). 

are thus U p /U c ~ 16 and 8.5 respectively. This means that a significant amount of energy (U p — U c ) is not directly 
used in the vaporization process, but is stored elsewhere on the timescale of r r . 

There are excitations of the system, such as phonons of different symmetry, but also potentially spin fluctuations, 
etc., that make up the difference between the condensation energy and the measured optical vaporization energy. Let 
us consider spin excitations first. The energy required to heat the entire spin bath from 4.5 K to T c for x — 0.1 is given 
by Um = Ji5i^ 0K Cm{T)cIT, where Cm{T) is the magnetic specific heat. Using the published value [26] of Cm{T) 
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for undoped La2Cu04 (Cm in doped Lai.gSro.iCu04 can only be smaller), we obtain Um — 80 mJ/mol (0.01 K/Cu). 
Clearly, the magnetic system alone is not capable of absorbing (U p — U c )/ks ~ 1-9 K/Cu, its heat capacity being too 

small by a factor of ~ 190. Making the same estimate for the lattice excitations, we obtain U L = J 4 \ K C p (T)dT ~ 77 
J/mol (9 K/Cu) for x = 0.1 (T c = 30 K) and 240 J/mol (28 K/Cu) for x = 0.15 (T c = 38 K), where C P (T) is the 
experimentally measured specific heat|25j. The phonon subsystem can thus easily absorb the excess supplied energy, 
with U L /U P ~ 4.5 for x = 0.1 (and 11.6 for x = 0.15). 

This observed discrepancy between measured U p and thermodynamically measured condensation energy, as well 
as the T- and ^-dependence of the superconducting state depletion process can be naturally explained within the 
Rothwarf- Taylor (RT) model in the bottleneck regime, where the pairing bosons are reaching quasi-equilibrium with 
the QPs [5] on the 1 ps timescale and share some of the energy supplied by the optical pulses. 

In the RT model, the pair-breaking time (which corresponds to the condensate vaporization time when T > Tt) is 
given by r^ 1 = rjy/l/A + (4A(0) + 2n(0))R/r] where the initial QP and boson densities are n(0) and N(0) respectively 
[7J[8]. For weak photoexcitation, when both n(0), N(0) < Ut, where the threshold density is defined as tit = n/R, 
r r is independent of J 7 , and 2r r = n~ l . For intense photoexcitation, when either n(0), iV(0) > Ut, T r strongly 
depends on T ' . A strong JF-dependence of r r is not observed in our data, which implies that LSCO is in the " weak" 
perturbation regime over our range of J 7 , and so n = l/(2r r ) ~ 0.5x 10 12 s~ . To estimate Ut, we take R = 0.1 cm 2 s _1 
measured by Gedik et al. in YBCO [10] and obtain a threshold density n-r = rj/R » 0.8 x 10 20 cm~ 3 « 0.8 x 10~ 2 /Cu. 
We can make an alternative microscopic estimate of n? using the formula for the bare recombination rate from ref. 
[2"7] (with phonons as the mediating bosons) R = ^3^ A N , where No is the density of states at the Fermi level, A 
is the superconducting gap, f2u is the characteristic phonon frequency and A the electron-phonon coupling constant 
(which is the same as appears in the McMillan formula for T c [27). Taking typical values for LSCO No — 5/eVCu, 
A as O.OleV , Hd ~ O.leV and the measured A = 0.9 [28], we obtain R = 0.7 x 10~ 8 cm 3 s _1 which gives a very similar 
threshold density as the phenomenological estimate n? — n/R — 1.5 x 10~ 2 /Cu. Note that both are just slightly 
lower than the estimated photoexcited QP density at threshold fluence which is n s p — — ~ 2.7 x 10^ 2 /Cu. We 
can conclude that the RT model involving phonons in the pair-breaking process gives a self-consistent quantitative 
description of the vaporization dynamics. 

Let us now see whether the relaxation processes on the sub- 1 ps timescale might somehow involve spin excitations. 
In this scenario, energy might be initally transferred from PE carriers to the spin subsystem on a timescale much 
shorter than 1 ps and QPs would then be excited from the condensate by absorbing energy from the hot bath of 
spin excitations. For energy relaxation only real (not virtual |20j ) processes are relevant and the relevant interaction 
between QPs and spin excitations is spin-orbit coupling. Such a scenario is consistent with our data, provided that 
the spin-orbit relaxation time ts-o is equal to, or shorter than the observed vaporization time of r r = 0.8 ± 0.15 
ps. To estimate the vaporization time for this case, we use the fact that spin-lattice relaxation is a process in which 
electron-phonon relaxation follows spin-orbit relaxation, and ts-l — T s-o + t e-p- So, for spin excitations to be 
involved in the pair breaking and QP relaxation process, ts-l needs to be of the order of 1 ps or less. Electron- 
paramagnetic resonance (EPR) measurements of Cu spin relaxation in Lai.gSro.iCuO^, give EPR lincwidths ranging 
from AH ~ 1 kilogauss at 30 K to AH = 3 kilogauss at 8 K. This corresponds to a lower limit of the relaxation time 
ts-l — 100 — 340 ps [23], which is much longer than observed. Assuming that the measured ts-l is correct, the pair- 
breaking thus cannot proceed via the spin excitations, because the relaxation process at 4.5 K would take over 340 ps, 
instead of ~ 0.8 ps. Thus spin excitations cannot be responsible for the destruction of the superconducting condensate 
by any currently known spin-orbit relaxation mechanism. This conclusion has important implications for the pairing 
mechanism in these compounds. The pair-breaking process discussed above is related to QP recombination (pairing) 
by time-reversal symmetry, and therefore both processes must involve the same mediating boson, i.e. phonons. We 
conclude that only phonon-mediated vaporization is consistent with the observed dynamics, effectively ruling out 
spin-mediated QP recombination and pairing in these materials. 

We wish to acknowledge valuable discussions and important comments from K. Alex Muller, N. Ashcroft, P. B. 
Allen, A. S. Alexandrov, D. Van der Marel, E. Maksimov, I. Bozovic and D. Newns. 



5 



I. SUPPLEMENTARY MATERIAL 



A. Calculation of the behavior of the photoinduced reflectance as a function of fluence T for a 

superconductor in the high density regime 

To accurately identify the point where the superconducting condensate is vaporized, we need to account for geo- 
metrical aspects due to the finite absorbtion length of pump and probe light as well as the transverse beam profiles. 

In the probe beam, the sample penetration depth has to be accounted for twice (upon entering and exiting the 
sample). The relative photoinduced change in reflectivity is then given by: 

^oc^ J exp ^-2^^cxp (-^jn qp (r,z)rdrdz (2) 

where X op is optical penetration depth, and 2p pr is probe beam diameter on the sample. We take the change of 
reflectivity to be linear with the density of photoexcited quasiparticles n qp , and n qp to be approximately linear with 
excitation density up to the threshold excitation fluence Tt where the superconducting condensate is evaporated. For 
T > Tt all quasiparticles are excited and uqp saturates at n s 



F(r,z) 

n„ 



T{r, z) < T T 
T(r, z) > T T 



(3) 



The light fluence penetrating into the sample is To = (1 — R)T pu where T pu represents the laser fluence on the surface 
of the sample. The laser fluence within sample is: 



T(r, z) = T exp 



exp 



P 2 

rpu 



(4) 



Assuming that in the normal state AR/R does not depend on T due to other processes, we calculate the integral 
in two parts: for T(r,z) < Tt and separately over the volume where Tt > Tt- With pump beam radius p pu and 

= \ -j — h -J- ) the complete integral is: 

Ptff \Plu P%r) 1 6 



To 
T, 



-n. 



3z 



exp 




e i- 



T(r,z) 



rdr 
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exp 



2z 



exp 



P 2 
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1 rdr 



dz 



(5) 



where 6 is Heaviside step function. The behavior of this integral depends on the ratio ^ . To get a dimensionless 
result we normalize the result with its saturated value: 



OO /"OO 



—jj- oc n s j j exp | 2- 
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For To < Tt the dependence is linear in T: 

AR _ 2T pj ff 
AR S ~ 3T T p 2 pr 



r 

~P 2 J 

rpr / 



2 \ 2 \ 

\ Ppr^op 

rardz = n s - J —. — 



(6) 



(7) 



while for To > Tt 
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6-3p 2 



(8) 
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FIG. 3: Amplitude as a function of fiuence. The p pu '■ p P r ratio is 1.5 for sample x — 0.1 and 4 for x = 0.15. The lines are best 
fits of the data to eq. [7] and [8] 

with p 2 = Pp U /pp r and / = As can be seen from eq. [7] one can determine the saturated value of excitation fiuence 
J~t simply by reading the value at which AR/R reaches 2 ( o 2 ^/3p 2 J , of the maximum (saturation) value of AR/R. 

To determine T vu in our measurements we used a pinhole of diameter 2r p i n and measured the power of the beam in 
front of the pinhole (Pj n ) and after the pinhole (Ptr)- With known repetition rate of the pulses (v rr ) we can calculate 
T pu and diameter of the beam p: 

T pu v rr f£ pm exp (-^27rrdr 

Ptr I Pin 7 (9) 

TpuVrr Jp 00 exp [ - ^ J 2-jrrdr 

and 

Pin = FpuV rr / exp ( - ^- ) 27rrdr = T heam -Kp 2 v rr (10) 



jo V P / 

To obtain the Tq we accounted for the light reflected from the cryostat window (8%) and the part reflected from the 
sample (R). 

The excitation energy density calculated from the absorbed energy fiuence is: 

U=^l (11) 

In our experiment we used the following parameters. The pump/probe beams radii were (42/28 (±15%)/xm) and 
(71/18 (±15%) pm) for the samples x = 0.1 and 0.15 respectively. Tq was determined with the relative accuracy of 
0.3. Penetration depths were calculated from the optical conductivity and the dielectric function data [30j and were 
found to be 147±15 nm for x = 0.1 and 156±15 nm for x — 0.15. The reflectivity is 0.135 for x — 0.1 and 0.130 for 
x = 0.15 [H]. 

The saturated fiuences obtained by fitting the eqs. 6 and 7 to the data (see Fig. [3]) are J 7 ? 1 =4.2±1.7/iJcm -2 
and J-"y 15 = 5.86 ± 2.3pJcm~ 2 for underdoped (x=0.1) and optimally doped (x=0.15) samples respectively. This is 
equivalent to the absorbed energy of 1.95±0.78/csK per Cu atom for x=0.1 and 2.56±1.0fcsK for x=0.15. 



B. The temperature dependence of the photoinduced change in reflectivity in the limit of condensate 

vaporization. 

In the following we describe the derivation of the expected temperature dependence of the photoinduced change in 
reflectivity at optical frequencies in the limit of excitation intensities higher than the condensate vaporization treshold. 
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Figure 2a) of the main text shows the time-variation of the photoinduced reflectivity trace recorded at 4.5 K on 
Lai.gSro.iCu04 (T C =30K) as a function of the photoexcitation fluence ranging from 1 to 200 ^J/cm 2 . Component 
A, which describes the dynamics of the superconducting state pair-breaking and recovery, exhibits clear saturation at 
high excitation intensities. Analysis of the T -dependence of the signal amplitude at t — 2 ps with the model described 
in the previous section gives the trcshold fluence Tt =4.2/iJ/cm 2 for x=0.1 and Tt =5.9/iJ/cm 2 for x=0.15. It follows 
from these data that at T > Tt the superconducting state is vaporized on the timescale of rs 1 ps after excitation 
with 50 fs optical pulses. 

The temperature dependence of the induced change in reflectivity at T > Tt is markedly different than the 
temperature dependence obtained in the low excitation regime, which is not surprising. Indeed, the amplitude of the 
induced change in reflectivity in the regime of condensate vaporization, A s , is given by 



AR 



R 



**>n -ff'S Rri Rs 



(12) 



where R n and R s are the reflectivities in the normal and in the superconducting states, respectivey, and R n ,R s 
R n — R s . At optical frequencies, the induced change in reflectivity is proportional to the induced change in the 
imaginary component of the refraction index and therefore proportional to the induced change in the real part of the 



optical conductivity ^ a " a cc giving 



A„ = 



AR 

~R 



(13) 



To determine the temperature dependence of A s we evaluate pr(T) where is given by the Mattis-Bardeen 
relation [3Tj 



S (M = 



(f(s)-f(e + tkj))g(s)d£ + 



(l-2f(e + huj))g(e)de 



(14) 



A—Hu 



Here / (e) is the Fermi-Dirac distribution function, hoj is the photon energy, A is the superconducting gap, and g (e) 
is 



9(e) 



e(e + Huj) + A 2 



Ve 2 -A 2 J(e + hw)-A 2 



(15) 



In the limit of A <C huo Eq.( 14 1 can be rewritten as 

ffiijj-A 



s i />nuj — i\ c\ />oo 

-i(fio;)~_ g( £ - hLj)d£+ —J^ f( £ ){g{s)-g{e-hw))de 



(16) 



The first term on the right hand side of Eq.( 16 1 can be evaluated exactly and is given bv[3T 



h= i 



Huj J \ 1 




(17) 



where E and K are the complete elliptic integrals. In the limit of photon energies being much higher than the gap, 
1^ < 1, it then follows that 



Huj \eA J 



(18) 



The second term on the right hand side of Eq.(16) can also be calculated exactly in the limit of T, A <§C Ku> and is 
found to be 



% /8¥Afc B T ( -A 
■ exp 



(19) 
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Clearly, I2 presents only a small correction which becomes noticable only in the close vicinity to T c 
is in the limit when fku 3> T, A given by 



^(T)oc^m 



Ahui 
iA(T) 



Therefore, A S (T) 



(20) 



Due to the fact that the data on the temperature dependence of the gap are fairly scarce, A(T) has been com- 
monly assumed to follow the mean-field (BCS-like) temperature dependence. Comparison of A S (T), where a BCS 
T-dependence of A(T) is assumed, to the experimentally measured | f ^|jr > ^r T shows, however, that A(T) has a 
substantially weaker T-dependence than BCS functional form. In fact, from the available data on the temperature 
dependence of A(T) in cuprate superconductors, which is obtained from the temperature dependence of the SIS 
tunneling junction caracteristics[32 , as well as from the intrinsic tunneling data[33 it follows that the gap follows 

A(T) = Ao ^1 — (T/T c ) 2 ^j temperature dependence over a wide temperature range. Indeed, using this functional 

form with huj = 1.5 eV, and 2Ao = 120 , a nearly perfect agreement between A S (T) and the experimentally mea- 
sured |^| JF>JF is found (Figure 2). We should note however, that unlike in the low excitation regime, where by 
fitting the temperature dependence of the amplitude of the photoinduced change in reflectivity the magnitude of the 
superconducting gap can be extracted, this is not the case in the high excitation regime. 



■ 20 uJ/cm 
o 40 uJ/cm 2 
A (Eq.(19)) 




FIG. 4: The temperature dependence of the experimentally measured |^p|jr >i7 r f° r J~ ~ 20 and 40 /iJ/cm 2 , compared to 
A S (T) given by Eq.(20l with A(T) = A (l - (T/T c ) 2 ) temperature dependence of the gap (solid line). 
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